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Abstract 

Acceleration of electrons by Langmuir plasma turbulence 
and by the Fermi mechanism are suggested to explain the 
two-stage acceleration processes observed in hard X-rays 
from solar flares. A model for particle acceleration in 
flares is presented. 
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Acceleration of Electrons in Solar Flares 


I . Introduction 

Recent observations of solar hard X-rays of energies greater 
than 10 keV by Frost and Dennis (1971) with an 0S0-5 instrument 
have shed new light on the particle acceleration processes in solar 
flares. The hard X-ray emission that they observed consists of an 
impulsive burst with photon energies up to 100 keV in the initial 
phase of the flare, and a second burst with much harder photon 
spectrum, extending up to 250 keV. Both types of burst showed 
nonthermal characteristics. The harder burst type maintained its 
power law spectrum for as long as 40 minutes. Frost and Dennis 
(1971) argue that these observations indicate the occurence of a 
two-stage acceleration process. In the first stage, which coincides 
with the flash phase of the Ho< flare, electrons are accelerated up 
to 100 keV on a time scale of only a few seconds. The second stage 
acceleration occurs almost immediately thereafter, and raises electrons 
to much higher energies. The rise time for the second stage burst is 
much longer than for the first, typically lasting a few minutes. 

Kane and Anderson (1970) in their observations of impulsive solar 
flare X-rays > 10 keV with instruments aboard 0G0-5 also find that 
the spectrum of the emitting electrons becomes much steeper beyond 
100 keV. However, they did not observe the second hard X-ray burst, 
probably due to instrumental limitations. 
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The idea of two stage acceleration has been inferred previous- 
ly by De Jager (1969) from solar radio burst observations. However, 
the mechanisms involved have not been definitely established. If 
the first stage is due to induced electric field (De Jager, 1969; 
Frost and Dennis, 1971) then it is difficult to explain why the 
induced electric field has a maximum value around 100 keV. A long 
standing problem is the establishment of an induced field in a high- 
ly ionized plasma as exists in the flare. Here we present an alter- 
native explanation: acceleration by Langmuir plasma turbulence that 
is generated by plasma instabilities. We show below that this 
mechanism is plausible in the initial phase of a solar flare. We 
believe that the electrons are continuously accelerated in the second 
stage and we show that the Fermi mechanism could be reponsible for 
the second stage acceleration. 

II . First Stage Acceleration 

In laboratory experiments, plasma turbulence is frequently 
observed when an instability sets in. Associated with such insta- 
bilities are the appearance of high energy particles, microwave 
emissions, and X-ray emissions (Smullin and Getty, 1962; Alexeff, 
Neidigh, Reed, Shipley, and Harris, 1963; Hamberger, Malein, Adlam, 
and Friedman, 1967; Babykin, Gavrin, Zavoisky, Nedoseyev, Rudakov, 
and Skoryupin, 1967). The situation is solar flares is rather 
similar. It seems resonable to investigate whether the-high 
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energy flare phenomena are due to plasma turbulence created by an 
instability. For a concrete example we consider the flare model 
in which annihilation of magnetic fields occurs at the neutral 
points or lines (Petschek, 1963) . In the vincinity of the neutral 
lines the large gradient in magnetic field B sets up a current 
density j as given the Maxwell equation 

V xB = j (1) 

V/ “ 


where c is the speed of light. Taking a plausible thickness 

S io^ cm for the diffuse region (Friedman and Hamberger , 1969) , 

10 —3 

we find for an electron density n 10 cm , AB/v/500 gauss, that 

the drift velocity of electrons, V^, is 


V 


d 


_j c A B 

n e 47Tn e ^ 


3 x 10 cm/sec 


( 2 ) 


where e is the electronic charge. This velocity is much greater 
than the thermal velocity of the electrons 

1 

2 



4 x 10 8 


cm/sec 


(3) 


6 

where k is Boltzman's constant, m the electron mass, and T ~10 K 

t? V 

the electron temperature in the lower corona. As the phase velocity 
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of the electron plasma oscillation, V lies in the range <£. 
V ph < OO , the condition of generating plasma waves via the Cerenkov 
emission process, V d >V th , * s well satisfied by the current-driven 
electrons moving with velocity V d « Note that electron plasma 
oscillations (also called Langmuir oscillation) with phase velo- 
city V ph ^ Vth are heavily Landau damped (Boyd and Sanderson, 
1969). As a result of this current-driven instability small scale 
Langmuir plasma turbulence is generated (Kadomtsev, 1965). Owing 
to the many degrees of freedom of a plasma, Langmuir plasma tur- 
bulence is also readily generated by shock waves , by high frequency 
radiations, and by many other plasma instabilities. 

From the above calculation we see that the Langmuir plasma 
turbulence can be generated by the flare instabiltiy. Once the 
plasma turbulence is generated, it is possible for the electrons 
to be accelerated statistically due to particle-wave interactions. 
The average increase of particle energy E with time t is given by 
the diffusion equation (Melrose, 1968): 


dE _ 1 J)_ 

dt e 2 a E 


(E 2 D) 


(4) 


The diffusion coefficient D is related to the energy distribution 
in phase space, w ( V ph)> of the Langmuir plasma turbulence by 
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where U) = (4TTne 2 /m ) 2 is the electron plasma frequency, V 

P© ^ 

is the electron velocity, and k = CJ pe /V p k i s "the wave number of 
the plasma wave. Notice that we have approximated the dispersion 
relation for electron plasma waves at finite temperature by UJ -C0 pe , 
where CO is the oscillation frequency of plasma waves. The lower 
limit in the integral is set at V^j because plasma waves with 
phase velocities less than V h are heavily Landau damped, and 
hence do not contribute to the interactions with particles. Now 
substituting (5) into (4) , we find for non-relativistic particles 

that 


dE 

cTF" 


e 2 CO 
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£ e W ( V 


m T 


ph 


= V) 


( 6 ) 


where m is the mass of the accelerated particle. The accelera- 
tion rate for non-relativistic particles depends therefore strong- 
ly on the turbulence wave spectrum. Since no waves with V ph < V th 
exist, we see from (6) that only particles with velocities v th 

are accelerated by Langmuir plasma turbulence. The energy of the 
plasma waves, W(V ph > , is related to the energy of the plasma tur- 
bulence, w(V ph ), by the relation 
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W < V ph> ® W < V ph ) 



(7) 


The spectrum of plasma energy density, w ( v ph ) > is not well known, 
and only a rough estimate is possible. Pikel'ner and Tsytovich 
(1968) have calculated the stationary spectrum of Langmuir plasma 
turbulence for the case T ~ T., where T. is the ion temperature, 
for excitation by many kinds of plasma instabilities. They show 
that owing to various non-linear processes, a stationay spectrum 
of Langmuir plasma turbulence can be produced. For example, in 
the current-driven instability discussed above, plasma waves with 
V ph V th are emitted initially. Then non-linear processes such 
as waves scattering on electrons, and on ions, and wave-wave inter- 
actions, transform the initially low phase velocity plasma waves to 
higher phase velocities. The shape of the resulting spectrum of 
the Langmuir turbulence depends strongly on the energy input, Q 
ergs/cm —sec, into the plasma oscillations. The bulk of the energy 
of the plasma turbulence will reside in a Maxwellian peak centered 
at the wave number k Q , if the energy input to the plasma oscillation 
Q is large. Now, 


k = k 
o * 


8 n T y 
e 

UJ Q 

pe w 


1 

2 V -2 


( 8 ) 


where k 

* 




(kT /4TTne^) 2 


is the Debye length, 
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\) ~ 50nT ” 5//2 is the electron collion frequency, and is a 

r e e 

numerical constant between 2.8 and 4. The condition for the 

formation of the Maxwellian bump is k Q k^. Otherwise the 

spectrum of plasma turbulence will be more or less flat within the 

range k . cz k* sS OJ/V . Taking T ~ 10 5 K, n ~ 10 10 cm -3 , 

and 3, we find the condition that the bump forms, k Q <<£ k*, 

13 3 

requires that Q » 4 x 10 erg/cm -sec. Even if we assume that 

32 

the total flare energy of ^ 10 ergs is dumped in the plasma 

25 3 

turbulence if a volume ^ 10 cm in one second, we will not obtain 
a large enough excitation Q to generate a bump in the Langmuir tur- 
bulence spectrum. In view of the great uncertainty of the spectrum 
of this turbulence, we shall assume that for the case of solar flares 
the resulting turbulence spectrum is constant, i.e., w = constant 
within the range k* ^ k ^ ^ pe / v th > and is zero outside. The 
acceleration rate (Equations (6) and (7)) is then 


dE 

” 


w 


m 00 


pe 


e 2 /yj 

m CG pe l 87C 


(9) 


where w = ~>|B 2 /8fC and is the fraction of the annihilated 

magnetic field energy that goes into plasma turbulence energy. 

10 -3 

Setting n ~ 10 cm , B ~ 250 gauss, we find, for electrons, m = m , 

© 



erg/sec . 


dE 

dt 


( 10 ) 
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In order that thermal electrons with E * 10 eV (T 10 K) be 

© 

accelerated, the acceleration rate must exceed the ionization 
loss rate given by 



( 11 ) 


where “ft = h/2ir and h is Planck’s constant. This requires that 

the transformation efficiency ^ ^ 6 x 10 -8 . Note that the thermal 

level of the Langmuir plasma turbulence w,./nkT ^ 1/(6 ft n An ) — 

3 x 10 -8 , giving ^ th = w th /(B 2 /8 7T) O: 2 x 10 -12 (for T 0 ^10 5 K, 

10 -3 

n 'v 10 cm , B '■v 250 gauss). The condition for weak plasma tur- 
bulence, w th «w«nkT e , is therefore well satisfied for 10 -4 
-9 

6 x 10 . Now the time for a thermal electron of energy ~ 10 eV 

to be accelerated to a final energy 100 keV is given by 



For 10" 4 ^ ^ 6 x 10 8 , the time needed for acceleration is 

—5 

10 ^ t ^ 0.3 second, which is consistent with the impulsive hard 

X-ray observations. The value of ^ in terms of the physical 
parameters of a specific instability is not known. However, Zaitsev 
(as quoted by Bhatia and Tandon ,1970) has estimated that in a 
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8 9 “3 4 

collisionless shock with n <*-' 10 - 5 x 10 cm } B ^ 300 - 10 gauss, 


-3 -4 

is around 10 - 10 


Our estimate of 


the maximum value of ^ 

is therefore not unresonable. 

We have shown that the Langmuir plasma turbulence is capable 
of accelerating electrons in a time scale compatible with the obser- 
vations. The cut-off in the electron spectrum is due to the cut- 
off in the spectrum of the plasma turbulence which has a minimum 

5 

wave number k* c* 0.6 (for T g 10 K) . This corresponds to V ^ ^ 
10^ cm/sec j or a particle energy of the order of 100 keV. Since 
w(Vpk 10^® cm/sec) ^ 0, no particles can be accelerated beyond 
about 100 keV. Note that the acceleration rate by Langmuir plasma 


turbulence is, in this case, inversely proportional to the particle 
mass. Hence protons are not as efficiently accelerated as electrons 
in the impulsive phase of solar flares. This agrees with Lin 
(1969) , who observed that in many cases only 10 - 100 keV electrons 
are observed in interplanetary space following solar flares, while 
no protons are associated with them. 


Ill . Second Stage Acceleration 

Frost and Dennis (1971) suggest that in the second stage the 
acceleration is associated with a shock front that also generates 
the type II radio burst. The type II burst at frequency of 80 MHz 

5 

for the March 30, 1969 event occurs at a plasma level •'v 4 . 5 x 10 km 

5 -3 

with a local electron density 8 x 10 cm (if we assume that the 
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density over the active region is 10 times greater than that of 
the normal corona) . Frost and Dennis (1971) suggest that the 
electrons that emitted hard X-rays were accelerated at about 18 
seconds earlier than the type II burst. If this is so, then, the 
level of the hard X-ray electron acceleration occurred at the height 
4.3 x 10 5 km, assuming that the shock front propagated at a speed 
of about 10 3 km/sec (Kundu, 1965). This height of the hard X-ray 
burst corresponds to a density of 10 cm" , which is much smaller 
than those usually deduced from solar X-ray observations (Takakura, 
1969; Holt and Cline, 1968). Thus, it remains to be shown that the 
shock front indeed accelerates the electrons that emit the hard 


X-rays . 

An important question is how the electrons could maintain their 
nonthermal power law for as long as 40 minutes, at least in the ob- 
served energy range from 40 keV to 250 keV while losing energy by 
both radiation and ionization. Let us consider an arbitrary an- 
isotropic electron distribution that is initially accelerated. 

Owing to electron-electron collisions with the ambient plasma, the 
high energy electrons are scattered in momentum as well as in energy 
space, and as a result the distribution tends toward isotropy and 
becomes Maxwellian at about the same time scale as given by (Rossi 
and Olbert, 1970) , namely: 


t ££ 
ee 



4 (21C) 2 e 3 n InA 


sec 


(13) 
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where T g is the final electron temperature, n is the electron den- 
sity, and InA is the usual Coulomb logarithm. Taking In A — 18 
in the ranges of density and temperature of interest to solar flares 

Q 1 O 4 7 

(n /s/ 10 - 10 cm , T ~ 10 - 6 x 10 K) , we find that t £ 7 

seconds. That is the electron distribution will become thermalized 
in a time scale much short than that indicated by the second hard 
X-ray burst. The long duration of the second hard X-ray burst with 
nonthermal characteristics indicates then that electrons are conti- 
nuously accelerated, and thereby replace the electrons which have 
been thermalized. It can also be argued that the electrons are 
initially accelerated by some mechanism, then trapped and stored 
in the corona. In that case, the hard X-rays would be produced 
through the gradual leaking of electrons from the storage region 
to lower levels of the atmosphere. However, it is not clear how 
a large number of high energy electrons could be confined by the 
relatively small magnetic fields in the corona. Here we suggest 
that in the second stage, particles are continuously accelerated by 
the statistical Fermi mechanism in collisions with the magnetic ir- 
regularities. These irregularities are continuously generated by 
the magnetic disturbances which are associated with the continuous 
energy input to the flare. 

The rate of energy increase for the statistical Fermi mechanism 


is given by 
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dE 



(14) 


where u is the velocity of the propagation of the magnetic inho- 
mogeneity, an 'T / is the mean collision time of the particles with 
the magnetic irregularities. In order for the acceleration to be 
effective the energy increase due to the Fermi mechanism must be 
greater than the ionization loss, which is the most important loss 
for electrons with energies less than 100 keV. From (11) and (14) , 
we find that the injection energy is 



The velocity of the propagation of the magnetic inhomogeneity is of 

1 s 

the order of the Alfven velocity, i.e., u as B/(47Tm n ) 2 -as 10 cm/sec 

10 -3 ” 

for B ~200 gauss, and n 10 cm (m p is the proton mass). Taking 
the injection energy to be 10 keV which is provided by the initial 
heating and the first stage acceleration, we find from (15) that 
the mean collision time of particles with the magnetic irregularities 

_5 

is T2; 2 x 10 second. This gives a characteristic scale of 

5 

field inhomogeneity of 1.5 x 10 cm. This length scale is of the 
same order of magnitude as that of the fine structure of inhomo- 
geneity observed in Ho(, flares (Suemeto and Hiei, 1959). Now, 
the time needed for acceleration from an initial energy E^ to a 


final energy E f is, from (14), 
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t o; 
acc 


T c' 
4 u" 


In 


E f 

E7 

i 


(16) 


g 

With 10 keV, E f ^ 250 keV, u ~ 10 cm/sec, the acceleration 

time is t Csr 2 seconds. Therefore, the Fermi mechanism is 
acc 

capable of accelerating electrons to high energies in a very short 
time interval. The final electron energy is limited by the brems- 
strahlung and the synchrotron radiation losses (Svestka, 1970). 

IV. Discussion 

In view of the above calculations, we propose the following 
model for acceleration of electrons in solar flares, with particular 


reference to the March 30, 1969 event (Frost and Dennis, 1971). 

The instability that triggers the release of the flare energy 
generates a Langnuir plasma turbulence. The plasma turbulence 
then accelerated electrons stochastically (cf.: Hall and Sturrock, 
1968; Tsytovich, 1966). The cut-off of the spectrum of plasma 
turbulence determines the cut-off of the acclerated electrons . The 
acceleration time is short enough to be consistent with the observa- 
tions, as we have seen. The electrons that escape outward then 
give rise to the type III radio burst while the electrons that 
move downward to the denser atmosphere emit the impulsive hard 
X-rays. After the initial instabilities and the collapse of the 
gases, hydromagnetic turbulence is generated. And we believe that 



-14- 


in this stage, stored flare energy is gradually released in a less 
violent manner than in initial instability of the flare. This 
release of the flare energy could occur through magnetic field 
annihilation (cf.: Petschek, 1963) . In the present model the 
period of bulk energy release is identified with the long time 
scale observed for second stage hard X-ray burst. In this stage, 
because of the continuous input of energy, magnetic disturbances 
are continuously generated. Hence the Fermi statistical mechanism 
continuously accelerates electrons and protons to high energies, 
giving rise to the proton events, the type II radio burst, and the 
gradual hard X-ray burst. At the same time, the ambient plasma 
is heated and gives rise to soft X-rays; thus a thermal plasma co- 
exists with a population of nonthermal high energy electrons. 



-15- 


References 

Alexeff , I., Neidigh, R. V., Reed, N. F. , Shipley, E. D. , and 
Harriss, E. A.: 1963, Phys . Rev. Letters, 10, 273. 

Babykin, M. B. , Gavrin, P. P., Zavoisky, E. K. , Nedoseyev, S. L. , 
Rudakov, L. I., and Skoryupin, V. A.: 1967, Soviet Phys.-JETP, 
25, 421. 

Bahatia, V. B. , and Tandon, J. N.:1970, Astrophys. Letters, 6, 113. 

Boyd, T. J. M. , and Sanderson, J. J.:1969, Plasma Dynamics, Barnes 
& Nobles, Inc., New York. 

De Jager, C. :1969, in Cospar Symposium on Solar Flares and Space 
Research , ed. by C. De Jager and Z. Svestka, North-Holland 
Publ. Co., pp 1 - 15. 

Friedman, M. , and Hamberger, S. M.:1969, Solar Phys., 8, 104. 

Frost, K. J., and Dennis, B. R.:1971, Astrophys. J., 165 , 655. 

Hall, D. E., and Sturrock, P. A.:1967, Phys. Fluids, 10, 2260. 

Hamberger, S. M. , Malein, A., Adlam, J. H. , and Friedman, M.:1967, 
Phys. Rev. Letters, 19, 350. 

Holt, S. S., and Cline, T. L.:1968, Astrophys. J., 154, 1027. 

Kadomtsev, B. B.:1965, Plasma Turbulence , Academic Press, New York. 

Kane, S. S., and Anderson, K. A.: 1970, Astrophys. J., 162, 1003. 

Kundu, M. R.:1965, Solar Radio Astronomy , Interscience Publ., New 
York. 

Lin, R. P. :1969, Acta Phys. Acd. Sci. Hungaricae, 29, Suppl. 2 
(Proc. 11th Conf. on Cosmic Rays, Budapest, 1969), 669. 



-16- 


Melrose, D. B.:1968, Astrophys. Space. Sci., 2, 171. 

Petschek, H. E.:1963, in Physics of Solar Flares , AAS-NASA symposium, 
NASA SP-50 , pp 425 - 439. 

Pikel'ner, S. B. , and Tsytovich, V. N..-1969, Soviet Phys.-JETP, 

28, 507. 

Rossi, B. , and Olbert, S.:1970, Introduction to the Physics of 
Space, McGraw-Hill Book Co. , New York. 

Smullin, L. D. , and Getty, W. D.:1962, Phys . Rev. Letters, 9, 3. 
Suemeto, Z., and Hiei, E.:1959, Publ. Astron. Soc. Japan, 11, 185. 
Svestka, Z.:1970, Solar Phys., 13, 471. 

Takakura, T. :1969, Solar Phys., (5, 133. 

Tsytovich, V. N.:1966, Soviet Phys . -USPEKHI , 9, 370. 



